On calculation of effective conductivity of inhomogeneous metals 
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Abstract 

In the framework of the perturbation theory an expression suitable for calcula- 
tion of the effective conductivity of 3-D inhomogeneous metals is derived. Formally, 
the final expression is an exact result, however, a function written as a perturbation 
series enters the answer. More accurately, when statistical properties of the given 
inhomogeneous medium are known, our result provides the regular algorithm for 
calculation of the effective conductivity up to an arbitrary term of the perturbation 
series. As examples, we examine (i) an isotropic metal whose local conductivity is a 
Gaussianly distributed random function, (ii) the effective conductivity of polycrys- 
talline metals. 
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g 1 Introduction 

1 — 'Rather often macroscopic properties of stochastically inhomogeneous media, which are 
i— h isotropic in average, can be described in the framework of different models of effective 
^ isotropic medium. Calculation of the effective conductivity (EC) is one of the well-known 
^ problems of theory of conductive media. It has been studied theoretically for many years 
(see, e.g., the review paper by A.G.Fokin However, in the general case up to now 
^ this problem has not been solved yet. 

Exact solutions for effective characteristics of stochastically inhomogeneous media 
^ ban be found very rarely. One of such examples is the calculation of EC of some two- 
c3 dimensional inhomogeneous media |2J]. The existence of this result is due to a specific sym- 
£3 metry transformation allowed by the equations of the problem. Recently in the frequency 
region of the local impedance (the Leontovich) boundary conditions applicability the ex- 
C3 act solution for the effective surface impedance of inhomogeneous metals was obtained 
q H |. Sometimes effective characteristics can be estimated qualitatively. For example, 
^ there are no regular methods allowing us to calculate the effective conductivity of three- 
dimensional (3-D) polycrystals accurately, but some geometrical reasoning suggested in 
^jjj Ref.0, allowed the authors to estimate EC of strongly anisotropic polycrystals. 

When the exact solution for an effective value ip e f of a characteristic if) cannot be 
found, perturbation theory calculations may be useful. In a stochastically inhomogeneous 
medium if> is a stochastic function of position r. A natural zero order approximation 
for ip e f is < if) > that is of ^(r) averaged over all possible realizations of the medium. 
However, when calculating < if> >, we do not take into account corrections due to the 
spatial fluctuations of i/j(t). 

From our point of view, the most accurate and physically meaningful method to take 
account of spatial fluctuations goes back to the work of I.M.Lifshitz and L.N.Rosenzweig 
. They proposed to start from differential equations of the problem coefficients of which 
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are random functions of position. Averaging these equations we derive equations for the 
averaged fields, which allow us to determine the effective characteristic. Usually this 
approach works when the inhomogeneity is small and perturbation theory is applicable. 
Previously it was used to calculate a lot of static and dynamic characteristics of inhomo- 
geneous media. As an example, we cite Refs. 0, ||, |9], [10], |TT|. In all these works only 
the first nonvanishing term taking account of spatial fluctuations was calculated. For a 
characteristic ip this term is of the order of < (ip— < ip >) 2 >. Usually, it does not depend 
on the correlations of the values of ip in different points of the medium. 

It is much more difficult to calculate the corrections of the higher orders in (ip(r)— < 
ip >)• The calculation involves step by step derivation of equations for each correction 
of the order n > 2. In addition, the high order corrections depend on the statistical 
properties of the medium. 

In this paper with the aid of the perturbation theory we derive an expression suitable 
for calculation of EC <r e / of 3-D inhomogeneous metals. The basic formula for <j e f is ob- 
tained in Section 2. Formally, in the framework of the perturbation theory applicability, 
the final expression, Eq.(15), is exact. However, a function written as a perturbation series 
enters the answer. So, more accurately, when statistical properties of the given inhomoge- 
neous medium are known, our result provides the regular algorithm for calculation of EC 
up to an arbitrary term of the perturbation series. We also write down the perturbation 
theory expression for EC of two-dimensional (2-D) conducting media. As an example, we 
examine an isotropic metal whose local conductivity is a Gaussianly distributed random 
function, and calculate the forth order correction to the effective conductivity. 

In Section 3 we consider polycrystalline metals. To verify our perturbation theory, in 
Subsection 3.1 we compare the perturbation theory calculation of a e f for 2-D polycrystals 
with the exact result of [0. In Subsection 3.2 EC of 3-D polycrystals is discussed. We 
calculate EC up to the third order term in the parameters of anisotropy and outline the 
region of the perturbation theory applicability. Also we examine the dependence of the 
third order term on the statistical properties of the polycrystal. 

2 Perturbation Theory for Effective Conductivity 
of 3-D Inhomogeneous Metals 

We consider a stochastically inhomogeneous metal that is isotropic and homogeneous in 
average. The elements of the local conductivity (LC) tensor a ik are stochastic functions 
of position r. By < ... > denote the ensemble average over all possible realizations of the 
medium. The averaged conductivity < cr ifc (r) >=< o > 5 ik , and 

a ik (r) =<a> (8 ik + A ik (r)), < A ik (r) >= 0. (1) 

The stochastic tensor A ik (r) describes the spatial fluctuations of LC. 
By definition, the effective conductivity a e f is specified by equation: 

<j(r)>=a e/ <E(r)>, (2) 

where < j > is the macroscopic direct current density and < E > is the uniform macro- 
scopic electric field; j(r) and E(r) are the local current density and the local electric field, 
respectively. 
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To calculate a e f, we write the electrostatics equations 

divj = 0, rotE = 0, (3) 

and the material equation that is the Ohm law. According to Eq.(l) 

ji(r) =<a> (5 lk + A ik (r))E k (r). (4) 

Because of the tensor A ik (r), Eq.(3) constitute a system of stochastic equations for the 
local field E^r). 

We write j(r) and E(r) as j =< j > +5j and E =< E > +5E (< 6j >=< 5~E >= 0); 
Sj and 5E are position dependent stochastic current and stochastic field, respectively. 
Substituting these expressions in Eq.(4) we obtain after averaging 

< ji >=< a> (<Ei> J, =< A ik (r)SE k (r) > . (5) 

Comparing Eq.(2) and Eq.(5), we see that the uniform vector J defines the contribution 
of the spatial fluctuations of LC to the value of EC. 
Subtracting Eq.(5) from Eq.(4) we have 

5j i =<a>(5E t + A ik <E k >+D t ), D { = A ik 5E k - J t . (6) 

The components of the vector D at least quadratic in powers of the elements of the tensor 
A ik . When calculating EC up to the first nonvanishing term taking account of spatial 
fluctuations, the vector D in Eq.(6) has to be omitted. However, just this vector defines 
corrections of higher orders. 

To calculate the vector J we use the perturbation theory. We present the stochastic 
field SEi(r) as a series in powers of the elements of the stochastic tensor A ik (r): 

oo 

^(r) = 5:^ B) (r). (7. a) 

n=l 

and seek the vector Jj as the series 

oo 

^ = E4" } , j} n) =<A jk (r)6E<r 1) (r)>. (7.6) 

n=2 

With regard to Eq.(6) from the electrostatic equations (3) it follows that the Fourier 
coefficient of the stochastic field 5Ei(r) is 5Ei(k) = —KiKj(Aj k {k) < E k > +Dj(k)), 
where Ki = ki/k. Consequently, the Fourier coefficients of 5E\ n \r) are 

5E ( i 1 \k) = -K l K j A jk (k)<E k > and ^ n) (k)| n>1 = -n^Df } (k), (8) 

where Dj n \k) is the Fourier coefficient of -Dj n ^(r) = Aj fc (r)5E^ n '^(r) — jj n ^. 

Equations (7) and (8) allow us to write the n-th term of the series (7.b) as a sum of 
the lower orders terms jj m ' ) (2 < m < n). Indeed, for n > 2 from Eq.(8) it follows that 

SE\ n \r) = - {// d'kA^A^AK - K-i)SEt'\K-i)^ - \j[ n) ) , (9) 

where = n k n \ 
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Next, we use Eq.(9) to write down the expression for 5E\ n ~ 1 \r). We Fourier analyze 
5E'j- n_1 ' ) (r) substitute the result in Eq.(9). Then the term proprtional to j/ n_1 ^ appears 
in the right-hand side of Eq. (9) and the Fourier coefficient SE[ n ^ (k) enters the integral 
term of this equation. 

The next steps are obvious. In successive order we decrease the superscripts in the 
expressions for the stochastic fields entering the integral term of Eq.(9). At the last 

step the coefficient 5E^(ki) defined by Eq.(8) enters the equation. 

If we use the expression for 5E\ n \r) obtained as the result of the aforementioned 
decreasing procedure when calculating J^ n+1 ^ and replace the superscript (n + 1) by m, 
we have 

m 

£4 m -°4° = -3t4 m) <£ fe >. (io) 

1=2 

In Eq.(10) wfj} = 5 ik , = - < A ifc (r) > /3 = and for n > 2 

»$? = ("D n i ( i ^nT /-/ d^-d-r^ /.../ d>k 1 ..A>k n - 1 $ 1 ...<fc^ 1 x 

< A iIl (r)A ilfa (r 1 )...A ifl _ lfc (r B _ 1 ) > e W-*) e tofr-**)... e *«-*n-*-*n-i) (n) 

When the medium is isotropic in average, the moments < A^ 1 (r)A J1 i 2 (ri)...Aj 2 / c (r n ) > 
depend only on the modulus of the differences between the vectors i"j entering the given 
average. Then the elements of the tensors do not depend on position. Since all 
isotropic uniform second rank tensors reduce themselves to the unit tensor 5^, the only 
possible form of all the tensors (n > 2) is 

W £> = w^ = lw%\ (12) 

where w^f is the trace of the matrix We also have = 1 and = 0. 
Taking account of Eq.(12) we write down Eqs.(10) for all 2 < m < oo: 

jf ) = -3w (2) < Ei >, if m = 2; (13.a) 

w {1) 4 2) + 4 3) = -3u> (3) <Ei>, if m = 3; (13.6) 
w WjW+ w WjW + jW = -3 w W < E i >, ifm = 4; (13.c) 

w (3) j(2) + ^(2) j(3) + w (l)jW + j(5) = _ 3w ( 5) <Ei> ^ iim = 5 {l3 d) 

and so on. We add together all these equations. Then according to Eq.(7.b) we obtain 

oo 

T 3 J t = -3(T 3 - 1) < Ei >, T 3 = ]T (14) 

m=0 

Consequently, Jj = 3(1/T 3 — 1) < ^ >. Finally, with regard to Eq.(2) and Eq.(5) we 
have 

a ef =<a>(3/T 3 -2). (15) 
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Together with Eq.(15) we would like to present the expression for EC of 2-D inhomo- 
geneous metals. Repeating the previous calculations with regard to the dimensionality of 
the problem, we obtain 

oo 

a e/ =Zf(2/T 2 -l), T 2 =Y,v {m) (16-fl) 



m=0 



where again = 1, d' 1 ' = and for n > 2 we have 



,(n) 

4(2?r 



x 



< A iJl (r)A il , a (r 1 )...A in _ li (r n _ 1 ) > e * 1 (r-r 1 ) e ,k 2 (r 1 -r 2 )_ e ,k„_ 1 (r„_ 2 -r„_ l); (16 . fo) 

r n and k„ are 2-D position and wave vectors, respectively. 

Equations (15) and (16) are the basic formulas we use when calculating EC within 
the framework of the perturbation theory. These formulas look rather simple. However, 
usually even for rather simple models of inhomogeneous conducting media we can calculate 
T 3 (or T 2 ) up to the desirable order in powers of the parameters of inhomogeneity only. 

It can appear that no significant simplification has been done. However, when the 
perturbation theory is applicable, Eq.(15) and Eq.(16) provide a regular method for cal- 
culation of the high order terms in the expression for EC. 

In the end of this Section as an example, we examine an isotropic metal with Gaussian 
inhomogeneities. This example is rather simple, since, first, for isotropic metals the 
expressions for are simplified, and, second, the form of n-point averages entering 
these expressions are known. 

We suppose that LC is o"jfc(r) =< o > (1 + A(r))5 ife , where A(r) is a strictly stationary, 
zero-mean Gaussian process |fL2"| . Then the two-point average is 

< A(r)A(n) >= A 2 W G (\r - n|), (17. a) 

where A 2 =< A 2 (r) >, and the Gaussian correlation function W G (r) and its Fourier 
coefficient W G {k) are 

3 

W G {r) = e~ r2 / a2 and W G (k) = - 3 e~ (ka)2/4 , (17.6) 

(2Vtt) 

a is the correlation radius. Next, the average of a product of an odd number of A(r) 
vanishes, and the average of an even number of A(r) is given by the sum of the products 
of the averages of pairs of the A(r)'s taken in all possible ways, irrespective of order, e.g. 

< A(r)A(n)A(r a )A(r 3 ) >= A\W G (\r - ri |)H/ G (|r 2 - r 3 |) + W G (\r - r 2 \)W G (\r 1 - r 3 |) 

+W G (|r-r 3 |)W G (|ri-r 2 |)). (17. c ) 

For this particular medium in the sum T 3 all the terms w {2n - 1] = 0. We use Eq.(ll) 
and Eq.(12) to calculate and w^. Then 



= ^-j d 3 kW G (k). 
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Being the value of Wc{r) when r = 0, the integral in right-hand side of Eq.(18) is equal 
to one. Thus, = A 2 /9. Note, Eq.(18) is valid not only in the case of Gaussian 
inhomogeneities, but for any random zero-mean function A(r). Consequently, always the 
first (quadratic in fluctuations of LC) correction to EC is independent of the form of the 
correlation function (see, for example, |fT3|). 

Our calculations showed that w^> = A 4 (7r — 20/9)/9. Then according to Eq.(14) and 
Eq.(15) up to the terms of the order of A 4 we have 

TP = 1 + I {A 2 + A\n - f)} , aif =< o > {l - 1[A 2 + A 4 (vr - 7 -)]} . (19) 

For Gaussianly distributed random function A(r) it is rather simple to calculate 
up to the higher order terms in A 2 . However, we failed trying to obtain a general formula 
for these terms. 



3 The Effective Conductivity of Polycrystals 

Polycrystals are widespread case of inhomogeneous media, where the inhomogeneity is 
due to different orientations of discrete single crystal grains. Anisotropic properties of 
each grain are described by tensor characteristics. If crystallographic axes of the grains 
are randomly rotated with respect to a fixed set of laboratory axes, the characteristics of 
the medium measured in the laboratory coordinate system are stochastic tensor functions 
of position. However, the invariants of the tensors are the same for all the grains and 
do not depend on position. Consequently, effective characteristics of polycrystals have 
to be expressed in terms of the invariants of the tensors that describe the phenomenon. 
When a polycrystal is isotropic in average (it is untextured), its effective characteristics 
are isotropic tensors. For example, EC of an untextured polycrystal is of/ = o e fb ik . 

3.1 Two-dimensional Polycrystals 

For 2-D polycrystal there is the exact solution for a e f found by Dykhne 0. We use the 
Dykhne formula to verify the calculations of Section 2, comparing the perturbation theory 
results with the expansion of the exact solution in powers of the parameter defining the 
anisotropy of LC. We do not find out a procedure allowing us to calculate an arbitrary 
term of the perturbation series. Therefore we only checked the first and the second terms 
taking account of spatial fluctuations. This exercise, being a test for our theory, outlined 
the way of calculation of high-order terms in the expression for EC of 3-D polycrystals. 

Let be the single crystal 2-D conductivity tensor, and let o\ and 02 be its principle 
values. We set 

*=k<*+<»), A = ^4, (20) 
2 (cri + er 2 ) 

where a is the mean conductivity, and the dimensionless parameter A (|A| < 1) defines 
the deviation of the principle conductivities from a: <J\ = a(l + A) and 02 = ^(1 — A). 
Being the first order invariant of the tensor a^, the mean conductivity is independent of 
position. Evidently, a =< o >. 

According to Dykhne, EC of a 2-D polycrystal does not depend on correlations between 
stochastic functions A i / C (r) at different positions: a e f = ^o\o<i. Note, o\<j<i being the 
determinant of the tensor a^, is its invariant. 
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We rewrite the Dykhne formula as 

i A 2 A 4 

a ef = a^{\ - A 2 ) =a(l — - — ...) (21) 

and compare the terms of the series (21) with the result obtained with the aid of Eqs.(16). 

To start the calculation we write the local conductivity tensor with respect to the 
laboratory coordinate system: o ik == &(5i k + A ik (r)), where 

Aj fc (r) = AZ ik (r); Z ik = (a a a k x - a i2 a k2 ), (22) 

and ai k — «jfc( r ) is 2-D rotation matrix that defines the orientation of the crystallographic 
axes of the grain containing the point r with respect to the laboratory axes. The elements 
of 2-D rotation matrix depend on one angle only. Let it be the angle ip. The value of ip 
is a stochastic function of position. 

With regard to Eqs.(16) when calculating EC of 2-D poly crystals up to the terms of 
the order of A 4 , we need to calculate the terms v ^ (n = 2, 3, 4) of the series (16. a) for T 2 . 
Consequently, first, we must define the second, the third and the fourth moments of the 
stochastic tensor function A ik (r). The general statistical properties of polycrystals were 
discussed in Ref. [[|. Here we discuss them briefly, paying attention to calculation of the 
high-order moments. 

The only property of the medium that affects the ensemble averages is the rotations 
of the crystallographic axes of the grains. Then ensemble average becomes the average 
over all possible rotations of the crystallites. We assume that in the ensemble the angles 
defining the rotations of different grains are statistically independent. When calculating 
the second order two-point correlator < A i fc(r)A; i (r 1 ) > entering the expression for v^ 2 \ 
there are two cases to consider: 

1. r and ri are in the same grain; 

2. r and ri are in different grains. 

We denote the probability of the case 1 by W^Qr — rx|). Then 1 — W 2 is the probability 
of the case 2. The fact that the probability W 2 depends on r and I - ! only through |r — ri| 
is the consequence of our assumption that the polycrystalline medium is isotropic and 
statistically homogeneous. 

Let ip n be the angle defining the orientation of the crystallographic axes in the grain 
containing the point r n . In the case 1 the rotations defined by the angles ip and ipi 
are identically equal. So, these both angles define one specific rotation. In this case 
the second order two- point correlator < Aj^i^A^rx) > reduces to the second order 
one-point correlator S ik ^u =< A^(r)A^(r) >. This is an isotropic second-rank tensor: 
Sik-u = S5 k i. As far as Ajj = 0, Ai k A k i = 2 A 2 is the only nonzero second order invariant 
of the tensor A im (r). Then S = A 2 . 

Because all the rotations are assumed to be independent, in the case 2 we have 
< Aj fc (r)A;j(rx) >=< A ik (r) >< A K (rx) >= 0. Consequently, the case 2 does not 
contribute to the value of the second order two-point correlator < A i( ! C (r)Aij(r 1 ) >. Thus 

< A ik (r)A H ( ri ) >= A 2 W 2 (\r - ri \)8 kl . (23) 

Next, it is evident that the three-point average < Aj fc (r)A; m (r 1 )A ni (r 2 ) > is not equal 
to zero only if all the three vectors r, r± and r 2 are in the same grain. Let us denote the 
probability of this event as W3 = Ws(\r — ri|, |r — r 2 |, — r 2 |). Then 

< A ifc (r)Aj m (ri)A ni (r 2 ) >= S ik .i m . ni W 3 (\r - n|, |r - r 2 |, |r x - r 2 |), (24) 
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where S ik] i m]ni =< A ifc (r) A im (r)A ni (r) > is the third order one-point average. It is clear 
that Sik-im-,ni is the isotropic forth- rank tensor. One can easily verify that in the 2-D case 
all the invariants of the tensor A ik Ai m A ni are equal to zero. Then, without calculations 
we can state that all the elements of the tensor S ik .i m . ni vanish. As a result, v ® = 0, and 
the third order correction to EC vanishes. 

Finally, it is clear that there are only two possibilities when the four-point average 

< Aj fc (r)A; m (ri)A n p(r 2 )A ? j(r 3 ) > is not equal to zero. Namely, 

1. Each of the pairs of the four vectors, but not all of them simultaneously, are in the 
same grain. 

2. All the four vectors are in the same grain. 

For the case 1 by W4([r a , rj,], [r c , rj) we denote the joint conditional probability for 
the vectors r a and r& to get in the same grain, and, simultaneously, for the vectors r c and 
Yd to get in some other grain. The probability W4([r a , r&], [r c , rj) excludes the possibility 
for all the four vectors to be in the same grain. Evidently, since all the rotations of the 
grains are independent, in the case 1 the four- point average is a product of two one-point 
averages of the second order. Next, by W 4 ([r, r 1; r 2 , r 3 ]) we denote the probability of the 
case 2. Then 

< Ai fe (r)Ai m (ri)A n p(r 2 )A 9i (r 3 ) >= W 4 ([r,ri], [r 2 ,r 3 ])S ikM S nPiqi +W A ([r,r 2 ], [ri,r 3 ])S ikjnp S lmtq 

+W 4 ([r, r 3 ], [ri, r 2 ])S ikiqi S lminp + W 4 {[r, r ± , r 2 , r 3 ])R ikMtnPiqi , (25.a) 

where 

S ikM = < A ik (r)Ai m (r) >, (25.6) 
Rik,im,n P ,qi = < A ife (r) A, m (r) (r) A n „(r) A qi (r) > . (25.c) 

Let us calculate the elements of the isotropic tensors Si k -i m and Ri k j m ,np,qi- Taking 
account of the symmetry properties with respect to interchanges of the indices, we see that 
the isotropic forth rank tensor Si k -i m must have the form Sik,i m — o,5i k 5i m +b(5u5 km +8i m 5 k i) . 
To determine the values of a and b we calculate the two independent invariants of the 
tensor A ik Ai m . They are the contractions A kk A u = and A ik A ki = 2A 2 . Then 

A 2 

Si k -im = —(—SikSim + Su5 km + 5 im 5 k i). (26. a) 

If we use the explicit form of the tensor A ik (see Eq.(22)), it is easy to see that in the 2-D 
case the sixth rank tensor 

R%k,lm,np,qi *S'jfc,gj*S'/m,np A 8 kq Si m ^ np . (26.6) 

Now we are ready to calculate the terms and of the sum T 2 . The same as 
when deriving Eq.(18), with the aid of Eq.(23) we obtain v ^ = A 2 /A. As usual, the first 
correction term is independent of the form of the correlation function W 2 {r). 

To calculate we substitute Eqs.(25) and Eqs.(26) in Eq.(16.b) for n = 4 and 
calculate the contraction entering the integrand of the left-hand side of this equation. 
Then 

- (4) = ^E* (9) - (27-) 
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where 

X(1) = / / / ^^^(r - " r 2)^(r 2 - r 3 )W 4 ([r, r,}, [r 2 , r 3 ]), (27.6) 

X W = -J— 4 J J J d 2 nd 2 r 2 d 2 r 3 J J^kA^-^^-^x 

[2(K^K^f - l]tf(n - r 2 )W 4 ([r,r 2 ], [n,r s ]), (27.c) 

and 

X^ = -^J J J d 2 n d 2 r 2 d 2 r 3 J Jd^h^-^^-^x 

(^^fSir, - r 2 ) [W 4 ([r, r 3 ], fr, r 2 ]) + W 4 ([r, n, r 2 , r 3 ])]. (27.d) 

Because of the presence of the 5- functions in the integrand of the expression for 
we have = W 4 ([r, r], [r, r]). Thus, r = = r 2 = r 3 , and, consequently, all the 

four vectors are in the same grain. However, according to the definition, the probabil- 
ity W 4 ([r a ,r;,], [r c , I'd]) is not equal to zero only if the vectors r a and r& are in the same 
grain, but the vectors r c and are in some other grain. Then X^ = 0. The same 
argumentation is valid when calculating X^ 2 \ Really, after integration over r 2 , the proba- 
bility W 4 ([r, ri], [ri, r 3 ]) enters the integrand. This probability is equal to zero too. Then 
X' 2 (I. 

As a result, only the third term X*- 3 -* of the sum (27.a) contributes to v^ 4 \ In the 
integrand of Eq.(27.d) the sum of probabilities W 4 ([r, r 3 ], [ri, r 2 ])+W 4 ([r, ri, r 2 , r 3 ]) defines 
the combined probability for the pairs of the vectors r, r 3 and r 1; r 2 be in the same 
grains not making difference between the cases when the first and the second pairs get 
in one or distinct grains. This combined probability can be expressed in terms of the 
probability W 2 . Namely, W 4 ([r, r 3 ], [n, r 2 ]) + W 4 ([r, n, r 2 , r 3 ]) = W 2 (|r-r 3 |)W 2 (|ri -r 2 |). 
When substituting this equality in Eq.(27.d) and performing the integrations, we obtain 
?;W = A 4 /8. 

Finally, with regard to Eqs.(18) up to the terms of the order of A 4 we have 

A 2 A 4 A 2 A 4 , x 

T 2 = l + — + — and a ef = a(l- — -—). (28) 

(compare with Eq.(21)). We see that within the accuracy of calculation EC is the same 
as the one defined by the Dykhne formula. 

The results of this subsection can be formulated as follows: we verified the perturbation 
theory of Section 2 and find out the forms of the correlators (see Eqs.(23), (24) and (25)) 
entering the integrands of the terms of the series T 2 . With regard to the dimensionality 
of the problem, the same form of the correlators must be used when calculating EC of 
3-D polycrystals. 

3.2 Three-dimensional Polycrystals 

Let a polycrystal be composed from 3-D single crystal grains. The local conductivity 
tensor is defined by Eq.(l). If the polycrystal is isotropic in average, < a >= a = 
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(a i + a 2 + cr 3 )/3, <7j (i = 1, 2, 3) are the principal values of the single crystal conductivity 
tensor. The principle values of the tensor A ik are 



Ai = ^ = (2^-^-03). A 3 = ( / 2 " 3 ~ (71 ~^ ) . (29.a) 

(0-1+O-2 + O-3) (0"1 + 0"2 + 0"3) (°~1 + °2 + 0"3) 

Since Ai + A 2 + A 3 = 0, the anisotropy of LC of a 3-D polycrystal is defined by two 
independent parameters. With respect to the laboratory coordinate system 

A ik = AiaaCKfci + A 2 a i2 a fc2 + A 3 a i3 a k3 , (29.6) 

where ai k is the 3-D rotation matrix, whose elements are stochastic functions of position. 

In what follows we calculate the third order fluctuation correction to EC of 3-D poly- 
crystals. In contrast to the 2-D case, this correction does not vanish. To perform the 
calculation up we need to define the terms and of the series T 3 (see Eq.(14) and 
Eq.(15)). 

We calculate the two-point average entering the expression for the same as in 
Subsection 3.1. With regard to the dimensionality of the problem < A^(r) Ajj(ri) >= 
D 2 W 2 (\r— ri\)Ski/3 where D 2 is the second-order invariant of the tensor A^: D 2 = A ik A ki . 
Then with the aid of Eq.(ll) and Eq.(12) we have = D 2 /27. Again the quadratic 
term of the series T 3 and, consequently, the quadratic correction to EC do not depend on 
the form of the correlation function. 

When calculating we use Eq.(24) to calculate the three- point average < A ik (r) A/ m (ri) A ni (r 2 ) >. 
This average does not vanish only if all the three vectors r, 1*1 and r 2 are in the same 
grain. As in the 2-D case, the elements of the isotropic forth rank tensor Si^im-ni have 
to be expressed in terms of the invariants of the tensor A ik Ai m A ni . However, in contrast 
to the 2-D case, only one of two independent invariants of this tensor is equal to zero. 
Namely, A ik A ki A u = 0, but A ik A kl A H = D 3 , where D 3 = Af + A 3 2 + A 3 3 . Then taking 
into account the symmetry properties of the tensor Sik-i m -,ni we have 

D 3 

Sik;lm;ni = — [~ 25fcn^m + 3(5fc;5 nm + 5 km 5 n l)\. (30. a) 

The probability W 3 for the three vectors r, r 1 and r 2 to be in the same grain can be 
expressed in terms of the probabilities W 2 . Indeed, let us choose one of the three vectors 
r, i"! and r 2 . Since all the choices are equiprobable, the probability to choose one definite 
vector is equal to 1/3. Let it be the vector r. Then the joint conditional probability 
of choosing the vector r, and finding vector ri in the same grain with the vector r, and 
finding vector r 2 in the same grain with the vector r is W 2 (\r — ri\)W 2 (\r — r 2 |)/3. Now 
it is easy to see that 

W 3 = ^(W 2 (|r- ri |)W 2 (|r-r 2 |) + W 2 (K^ 

(30.6) 

When we take into account Eqs.(30) and substitute Eq.(24) in Eq.(12), after integra- 
tion over position vectors and over the angles in the double integral over the wave vectors 
ki and k 2 , we obtain 

w^ = -^\^--j\, (31.a) 
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j = 2tt 2 / / dk 1 dk 2 w 2 (k 1 )w 2 (k 2 ) i - , 1 ln K fc i + - fc 2 ))| 



'o Jo [ 2kiq 2 (kf + fc| 

Then up to the third order term in the anisotropy we have 
T 3 = 1 + w (2) + w {3) , a ef = a 



D, + £,n 1 _ J 

9 135 V 3 



(31.6) 



(32) 



From Eq.(32) it is clearly seen that EC of a 3-D polycrystal depends on the statistical 
properties of the medium. Of course, this is the well-known result. However, now we can 
estimate the influence of the form of the correlation function on the value of a e f. To do 
this we evaluated J for several choices of W 2 (r). The obtained results were 

1. J = 0.028, when W 2 = exp(-r 2 /a 2 ); 

2. J = 0.136, when W 2 = exp(— r/a); 

3. J = 0.052, when W 2 = 1/(1 + r 2 /a 2 ). 

For all the examined forms of the correlation function the value of J is very small. Con- 
sequently, although the value of the third order term in the expression for EC depends 
on the form of the correlation function, this dependence is very weak. 

Concluding this subsection, let us estimate when the low order perturbation theory 
calculations are sufficient. In Ref.[I3| it was shown that the effective conductivity is 
restricted by inequalities < j 2 > /j<3" _1 j < e / < eae/ < e 2 >, where <j -1 is the inverse 
of the local conductivity tensor <j; j and e are arbitrary test vectors of dimension of the 
current density and the electric field, respectively. If as the test vectors we choose a direct 
current j and a uniform electric field e , we obtain 

< a ef < < a > . 33 

< p > 1 

In polycrystals < cr > is equal to the mean conductivity {a\ + cr 2 + cr 3)/3 and < p > is 
equal to the mean resistivity (erf 1 + ovT 1 + Og )/3. 

For simplicity, let us examine a uniaxial metal: <j\ = a 2 . According to Eq.(29.a) in 
this case A x = A 2 = A = {a x - <t 3 )/(2<ti + a 3 ), and A 3 = -2A. Then -1 < A < 1/2, 
when < ai/o-3 < oo. In Fig.l following Eq.(33) we plot the lower and the upper bound 
of the ratio cr e f/a as function of A as well as the values of <J e f/a calculated with the aid 
of Eq.(32). If —.2 < A < .2 (or 4/7 < 01/03 < 2), the difference between the lower bound 
for EC and its upper bound is less than 10 percents. It is reasonable to assume that the 
perturbation theory results are applicable not only if 1 01/03 — 1| <S 1, but at least for all 
01/03 from the aforementioned interval. The Table 1 based on date from [15j shows that 
rather often the perturbation theory calculations are sufficient when estimating EC. 

We would like to note that when comparing the calculated value of a e f with exper- 
imental results, one must have in mind that in our calculations only the inhomogeneity 
due to different orientations of polycrystalline grains was taken into account. Of course, 
in real polycrystals there are other sources of inhomogeneity too. For example, we do not 
take into account the real structure of the boundaries of the grains. This simplification is 
justified when the grains are sufficiently large and the properties of the boundaries of the 
grains do not affect the result significantly. 
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List of Figures 

Fig.l. The lower bound (dashed line) and the value of o e $l~v (full line) for uniaxial 
polycrystals as functions of the parameter A. Always the ratio o e jja < 1. The value of 
cr e f/a is calculated with the accuracy up to A 3 . 
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Table 1. 



Met all 


Components of 
conductivity 
tenzor xlO 6 
(Ohm x cm) -1 


A 


<7 e# /t7 


Beryllium 


<J 1 = <7 2 =0.319 
<J 3 =0.279 


0.044 


0.999 


Bismuth 


<J 1 = <7 2 =0.0092 
<7 3 =0.0072 


0.078 


0.996 


Galium 


ST — A AC70 
U ^ — U.UJ / O 

a 2 =0.1274 
<7 3 =0.018 


Aj =-0.1467 
A 2 =0.8809 


0.85 


Cadmium 


O x = (J 2 =0.1475 
<7 3 =0.1196 


0.067 


0.997 


Magnesium 


<J 1 = <J 2 =0.221 
<7 3 =0.2646 


0.062 


0.997 


Tin 


<J 1 = <7 2 =0.101 
<7 3 =0.07 


0.114 


0.991 


IV ft 1 V\ 7 

iviercury 

(-45°C) 


(7 l — (7 2 — 0.043 
<7 3 =0.056 




QQS 
\J.yyo 


Stibium 


CJj — C> 2 — U.U24 

<7 3 =0.028 


-0 061 

\J .\J\J J- 


998 


Zink 


CTj = G 2 =0.172 
<J 3 =0.165 


0.012 


0.999 



